Coupling strength estimation for spin chains despite restricted access 
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Quantum control requires full knowledge of the system many-body Hamiltonian. In many cases this infor- 
mation is not directly available due to restricted access to the system. Here we show how to indirectly estimate 
all the coupling strengths in a spin chain by measuring one spin at the end of the chain. We also discuss the 
efficiency of this 'quantum inverse problem' and give a numerical example. 
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I. INTRODUCTION 

The great progress in experimental techniques for manip- 
ulating microscopic objects has brought even quantum me- 
chanical systems under control. It has been one of the strong 
driving forces to promote the recent intensive theoretical study 
of quantum information science. What we need towards the 
realization of quantum information processing and quantum 
simulations is control over a quantum system that resides in a 
high-dimensional Hilbert space. Then, in order to achieve the 
physical controllability, we typically decompose the Hilbert 
space into a tensor product of many low (usually two-) di- 
mensional spaces and consider addressing each of them indi- 
vidually as well as coupling any two members. 

Although the addressability of the whole system is thus 
the key, paradoxically it is the origin of the biggest problem 
in quantum control. That is, that we can access the system 
means that the surrounding environment can also interact with 
it, causing unwanted errors. In order to circumvent or mini- 
mize the effect of noise, various methods have been proposed, 
e.g., quantum error correcting codes, decoherence-free sub- 
spaces, and topological quantum computing. Here, we focus 
on an alternative approach, i.e., operating only a (small) sub- 
set of the system while isolating the rest from its environment. 
If we have all information on the Hamiltonian describing the 
system, its controllability can be analyzed by the theory of 
quantum control (see QJ for an introduction). Recently, it has 
also been shown that entire spin chains can be fully controlled 
by operating on one end only, provided all parameters of the 
Hamiltonian are known |2, 3]. 

However, if we could completely isolate part of the network 
from any interactions with the outer world, then this isolation 
would prevent us from obtaining the required information on 
the shielded subsystem. While the type of interaction that gov- 
erns the dynamics inside the network would be known due to 
its intended design, the precise values of various parameters 
that characterize the full Hamiltonian might only be in a cer- 
tain range and not accurate enough to enable us to have a de- 
sired controllability. Hence, a natural question arrises: is it 
possible to estimate the necessary parameters of the Hamilto- 
nian only by accessing a subset of the whole network? We an- 
swer this question positively and show how to do it for chains 
of spin- 1/2 particles whose interaction is of Heisenberg type. 
The main task is thus to estimate the coupling strengths be- 
tween 'untouchable' spins. This is also useful for quantum 
state transfer in spin chains, where some schemes require a 
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Figure 1: (Color online) All coupling strengths (black lines) of a 
chain of spins can be estimated indirectly by quantum state tomogra- 
phy on one end (dashed area). 



good knowledge of the system Hamiltonian f?) . 

The question of system identification of a quantum device 
by Fourier analyis has recently been studied in Jsj], where the 
identification of two-level subspaces was emphasized. Here 
we give a solution to system identification for the N-Ievel case 
under strong local constraints. Our work can be seen as an ex- 
ample of inverse problems that are very important in a num- 
ber of areas of science and engineering. There are plethora of 
situations where indirect probing is the only way to acquire 
desired information, from medical ultrasonography to seismic 
reflection in geophysics. A closely related problem to the one 
in this paper is the estimation of spring constants in classical 
harmonic oscillator chains |6]. Our contribution is to provide 
the quantum mechanical counterpart of this problem and to 
discuss the efficiency of the estimation procedure. Roughly 
speaking, we obtain information on the couplings by inserting 
one excitation into the chain and observing its return proba- 
bility. Since the excitation travels back and forth the chain, it 
obtains knowledge on all the couplings. This is proved by the 
main theorem in Section |II1 It turns out that in order to ob- 
tain good knowledge on the coupling strengths, the excitation 
needs to 'see' each single link of the chain 2N times, where 
N is the length of the chain. This claim will be made rigorous 
by using the uncertainty principle of the Fourier transform in 
Sectionlml 



II. SETUP AND MAIN RESULT 

We consider a chain of N spins coupled by an anisotropic 
Heisenberg Hamiltonian, i.e. 

iV-l 
n=l 

with unknown couplings (5„ and anisotropy A. Additionally, 
we assume that we know the signs of 8\ and ^2- ^ conserves 
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the number of excitations, i.e. 



.v£V 

This allows us to speak about local excitations in the system. 
In particular, we will see later that for almost all cases it suf- 
fices to restrict the parameter estimation to the single excita- 
tion sector of the system. As usual, we let |n) denote single 
excitation states with spin n in the state 1 1) and all other spins 
in the state |0). The state with all spins in |0) will be denoted 
as |0). The unknown part of the Hamiltonian is the interac- 
tion strength (5„ between spin each neighboring pair of spins 
(n, n + 1). The purpose of the following will be to estimate 
these coupling strengths. 

Theorem: Assume that the Hamiltonian ofEq. ((7) has a non- 
degenerate spectrum in the first excitation sector Then 
the coupling constants 5n can be obtained by acting on 
on the first spin only. 

We remark that although the non-degenerate spectrum is the 
generic case, in practice almost-degenerate eigenvalues can be 
problematic to the efficiency of the method suggested below. 
This will be discussed in more detail in Section |III] The fol- 
lowing lemma is similar to the inverse problems for classical 
oscillator chains considered in [6]: 

Lemma: Assume that all eigenvalues Ej (j = 1, . . . , N) in 
the first excitation sector of the Hamiltonian ofEq. (|7} 
are non-degenerate and known. Assume that for all or- 
thonormal eigenstates \Ej) in the first excitation sector 
the coefficients are known. Then the coupling 

constants Sn are known. 

While the requirements of the lemma may sound unrealistic at 
the first sight, we will see towards the end of this section that 
they are provided by a simple Fourier analysis of the return 
probability of a single excitation. 

Proof of the Lemma: Our first observation is that (setting 
6o = Sn = for the boundary terms) 



H\n) = S„-i\n-l) +Sn\n- 



(2) 



with 



= A{G-2Sn-2Sn-i) 

N-l 
G = 5ra. 

The first equation we use is 

D„ = {n\H\n)=J2E,\{n\Ef, 



(3) 



In particular this implies that _Di is known by the requirements 
of the lemma. Then, for all j we have 



Taking the inner product with (n| and using Eq. (|2]i we obtain 

{Ej - Dr,){n\Ej)-Sn-i{n-l\Ej) ^ 5r,{n+l\E,). (5) 

For n — 1 this reads 

{E,-D,){1\E,)=6^{2\E,). 

Since the l.h.s. is known for all j, the expansion of |2) in the 
basis \Ej) is known up to the unknown constant Si . Through 
normalization of |2) we then obtain \5i\ and since the sign 
of 6i is known by assumption, we obtain 6i. Also, we now 
know {2\Ej) for all j. The next equation is obtained by setting 
n = 2 in (|5]l as 



{E, - D2) {2\Ej) - 5i{l\E,) = S2{3\E^ 



31- 



The only unknown on the l.h.s. is D2 which is obtained 
by Eq. (|3]l. Using again the normalization we obtain ^2 and 
(3|£'j). We could continue this procedure, but the normaliza- 
tion would not provide us with the signs of the (5„ . We know 
more: from Di — D2 = — 4A(52 we obtain A and from Di we 
can then get X]m=i ^m- Then, Eq. (O gives us and there- 
fore including its sign. This method is then easily general- 
ized by using Eq. (|5]l for n = 4 to obtain {3\Ej), and Eq. ([3j 
for 64 , and so on. ■ 

Now we describe how the requirements of the lemma can be 
measured by controlling the first spin only. First, we initialize 
the system in ( 1 0) + 1 1) ) . We remark that this can be done 
by acting on the first spin only, cf. [f]). We then perform 
quantum state tomography ISJ on the same spin at a later time. 
The evolved state is 



-l=U{t) (|0) + |1)) = i= |^exp{-zGO |0) + J2 /ni Win) 

with fni given by (n|[/(i)|l). The reduced density matrix of 
the first spin of the chain is given by 



Pi 



2-l/nl' exp{iGt}/i*i 
exp{-*G<}/n l/nl' 



/2. 



E,\E,)=H\E,} 



(4) 



Thus, by performing quantum state tomography on spin 1 we 
can sample the following matrix element of the time evolution 
operator, 

exp{-GO(l|C/(t)|l) = E e^P{-* (^^- + ^) ^> • 

j 

(6) 

Since the eigenvectors of H are determined only up to an ir- 
relevant phase, we can assume, without loss of generality, that 

{l\E,)>0{j^l,...,N). 

Hence through Fourier transforming Eq. (|6|l, the coefficients 
(1 l-Ej) are known (as long as the spectrum is non-degenerate). 
The eigenvalues Ej are known up to the constant G. As usual, 
this global shift of eigenvalues does not change the physics - 
Eq. Q shows that G cancels out when determining the cou- 
pling strengths. 
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Figure 2: (Color Online) Simulated measurement data and its Fourier 
transform. The green stars show the position of the exact eigenfre- 
quencies. The coupling strengths have been chosen randomly in the 
interval 1 ± 0.05 for a chain of A'^ = 20. The simulated measurement 
points allowed the estimation with a standard deviation of 0.01 with 
respect to the real couplings. The Fourier transform was computed 
using standard FFT algorithms and a Hann window. Only the 10 
peaks with positive Ej are shown, the others are symmetric around 
the y-axis. 

III. EFFICIENCY 

The efficiency of the coupling estimation can be studied us- 
ing standard properties of the Fourier transform (see |9] for an 
introduction). The functions is sampled for each 

time t by state preparation, system evolution, and quantum 
state tomography. Therefore an important cost parameter is 
the total number of measured points, being proportional to 
the sampling frequency. The minimal sampling frequency is 
given by the celebrated Nyquist-Shannon sampling theorem 
as = Ema\, where i?max is the maximal eigenvalue 

of H in the first excitation sector. Due to decoherence and 
dissipation, the other important parameter is the total time in- 
terval T over which the functions need to be sampled to obtain 
a good resolution. This is given by the (classical) uncertainty 
principle that states that the frequency resolution is propor- 
tional to 1/T. Hence the minimal time interval over which 
we should sample scales as T^jjj = 1/ (Ai?)jj^jjj , where 
(Ai?)jjjjjj is the minimal energy distance of the eigenvalues of 
the Hamiltonian. Another important parameter is the height of 



the peaks in the Fourier transform given by |(l|£;j)r ■ These 
should be high enough to resolve, which means that all en- 
ergy eigenstates need to be well delocalized. If there is too 
much disorder, localization will take place (see for example 
JTqI ), and couplings far away from the controlled region can 
no longer be probed (in turn, this suggests a way of obtain- 
ing information on localization lengths indirectly - see con- 
clusion). When localization is negligible, the numerical algo- 
rithm to obtain the coupling strengths from the Fourier trans- 
form is very stable |6]. The reason is that the couplings are 
obtained from a linear system of equations, so errors in the 
quantum state tomography or effects of noise degrade the es- 
timation linearly. In our numerical analysis we found a good 
agreement with the real couplings for systems with small ran- 
domness (Fig. |2|i. Let us also look at the scaling of the prob- 
lem with the number of spins. Typically the dispersion rela- 
tionship in ID systems of length N is cos which means 
that the minimal energy difference scales as « -^^ and the 
total time interval should be chosen as N'^ . This agrees well 
with our numerical results (tested up to = 100). For each 
sampling point a quantum state tomography of a signal of an 
average height of « 1/A^ needs to be performed. Since the er- 
ror of tomography scales inverse proportionally to the square 
root of the number of measurements, roughly N'^ measure- 
ments are required for each tomography. 

IV. CONCLUSION AND OUTLOOK 

In conclusion, we have found that a vast class of spin chain 
Hamiltonians can be estimated by some restricted operation 
at the chain end. There are obviously many variants of the 
inverse problem that we have introduced above. For example, 
it is straightforward to see that the anisotropy parameter A 
can be estimated even if it is site dependent (i.e., A„) if one 
assumes that the signs of all (5„ are known. We have empha- 
sized here only one example as we think it is the most relevant 
one for the applications |3, 4] and stands as a paradigm for 
the setup we have introduced: quantum system identification 
of a black box by restricted access. There are many interest- 
ing more general and fundamental questions that our study 
suggests. For instance, is it possible to count the number of 
qubits by restricted access [11]? What can be said about more 
complicated networks and higher dimensional lattices? If the 
black box contains some dissipation and decoherence, what 
can we learn about it through restricted access? What can be 
learned about the localization length? 
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